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We use the replica trick approach to calculate, in the large-iV limit, the entanglement 
c-function in two-dimensional quantum chromodynamics with the gauge group U(N ) and 
quark in the fundamental representation (’t Hooft model). We show explicitly that the 
result is equal to a sum (with certain numerical coefficients) of entanglement c-functions 
of free massive scalars, which represent mesons in spectrum of the ’t Hooft model. Inter¬ 
estingly, each meson contributes as a point-like particle. We provide an explanation for 
why this is the case. 
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1. Introduction 


The subject of entanglement entropy in quantum field theory received a considerable 
attention over the past several years. A general physical setup is the following. One takes 
a system in the vacuum state 10), and specifies a sub-region A in the coordinate space. 
This way the whole system is split into sub-region A and its complement A . 1 The state of 
the subsystem A is described by a density matrix p. obtained by tracing out the degrees 
of freedom of the subsystem A. 

p = iv 3 |0){0|. (LI) 

The problem is to find corresponding entanglement entropy 


S — —Tr p log p , 


( 1 . 2 ) 


where trace is taken over the states of A. 

One way to calculate entanglement entropy (1.2) is a replica trick [1,2,3]. The simplest 
system, which this method can be applied to, is a single interval in 2d conformal field theory. 
If L is a length of the interval, and c is a central charge of the CFT, then entanglement 
entropy is given by [1,2,3] 

Sc FT — - log — , (1.3) 

3 e 

where e is a UV cutoff. The entanglement c-function for single interval in 2d field theory 
is defined as 


r dS 

Ce = l JL’ 


(1.4) 


and due to (1.3), for CFT it is equal to 


c 



(1.5) 


Moving away from the conformal regime quickly complicates calculation of the entan¬ 
glement entropy. The question which one can ask is what is the entanglement entropy of 
a single interval in the theory of free 2d scalar, or fermion, with a mass m. In both cases 
entanglement entropy is expressed in terms of solution of the Painleve V equation. This 
solution is known numerically for general mL, and is known analytically only in asymptotic 
UV and IR regions [4,5,6]. 

1 Below in this section we discuss subtlety of such splitting in systems with gauge interaction. 
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It becomes more challenging to find entanglement entropy in interacting quantum field 
theory. One of the computational tools has been suggested by Ryu and Takayanagi [7,8], 
and it gives a prescription to find entanglement entropy in strongly interacting quantum 
field theory, which has a known gravity dual. The other, purely field theoretical, method 
involves perturbative computation of entanglement entropy, see, e.g., [9,10] and references 
therein, for recent developments. 

In this paper we use the replica trick approach to calculate entanglement c-function in 
two-dimensional quantum chromodynamics with a single-flavored quark in the fundamental 
representation of the gauge group U(N). This system is known as’t Hooft model [11]. It 
is solvable in the large-TV limit. The spectrum of the’t Hooft model is given by a tower of 
mesons, with masses of sufficiently heavy mesons approximated by the formula 

m^~7r 2 A/c, (1.6) 

where A is the ’t Hooft coupling [11,12], see [13] for a review. The mesons are non¬ 
interacting in the large -N limit. 

Entanglement entropy in pure U(N ) gauge theory in two dimensions has been cal¬ 
culated and discussed, e.g., in [14,15,16]. There is a subtlety in defining entanglement 
entropy for systems with gauge interaction. One way to see the problem is to look at 
gauge-invariant operators of a finite size. The question is, how to treat the gauge-invariant 
operators which are partly located in A. and partly in A. Following [14,15,16] we as¬ 
sume that entanglement entropy is well-defined in 2d field theory with gauge interaction, 
and replica trick calculation is applicable, and gives a sensible answer. (See also [17] and 
references therein for discussion of entanglement entropy in gauge theories.) 

Each meson in the spectrum of the ’t Hooft model can be viewed as a string of a 
given tension. Length of such a string increases as mass of the meson grows. The subtlety 
in defining entanglement entropy arises, e.g., when we look at string which has one end 
inside the interval A, and the other end outside of the interval. At the same time, due 
to having a finite length, one would expect that mesons contribute to the entanglement 
entropy as finite-size objects. However, our replica trick calculation shows that each meson 
contributes to the entanglement c-function as a point-like particle. We explain that this is 
a consequence of a consistent definition of [16] of entanglement entropy in gauge theories. 

In two dimensions entanglement c-function (1.4) is a monotonically decreasing function 
of the (dimensionless) length of the interval [18]. Therefore it can be compared with the 
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Zamolodchikov c-function. The latter is defined via the two-point function of the stress- 
energy tensor [19,20]. Schematically, it is given by 

C(\ z \) — z 4 (T(z)T(0)} . (1.7) 

In’t Hooft model we have qualitative behavior of Zamolodchikov c-function given by 2 

C(V\\z\) = N F{V\\z\), (1.8) 

where F is a smooth function of order one, monotonically decreasing between one in the 
ultra-violet and zero in the infra-red. At a fixed point Zamolodchikov c-function is equal 
to central charge of the corresponding CFT. Similarly, entanglement c-function of CFT is 
proportional to the central charge (1.5). 

Due to these similarities between entanglement and Zamolodchikov c-functions, we 
expect the entanglement c-function for the ’t Hooft model to behave in a way, similar to 
(1.8). However, our calculation demonstrates a different behavior. In agreement with [18] 
we have obtained a monotonically decreasing entanglement c-function. We have derived 
that enanaglement c-functions is Ce(VXL) = 0( 1) for XL 2 3> 1 /N. unlike the Zamolod¬ 
chikov c-function (1.8). When XL 2 = 1/N, we obtained that the entanglement c-function 
approaches the value of O(N). In this regime interaction of quarks is parametrically sup¬ 
pressed, and perturbative calculation of entanglement entropy should be done. 3 We leave 
this for future work. 

As a warm-up, we start by analyzing a simpler theory with four-fermionic interaction, 
the Thirring model. This model is well known to be equivalent to a free massless scalar. 
This can be immediately demonstrated by bosonization of fermions. Therefore the single¬ 
interval entanglement entropy in Thirring model is given by the CFT formula (1.3) with 
central charge c = 1. We choose to postpone bosonization, and apply the replica trick 
approach to calculation of entanglement entropy in the Thirring model. We explicitly 
show that CFT result (1.3) for free complex-valued fermion (c = 1) is not changed by 
Thirring four-fermionic interaction. 

2 See App. B of [21] for results on two-point functions of general bi-fermionic operators in the 
’t Hooft model. 

3 See, e.g., [22,23] for examples of perturbative replica trick calculation of entanglement entropy 
in 2d quantum field theory. 
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Before proceeding to the ’t Hooft model, we study a simpler gauge theory in two 
dimensions, Schwinger model of a fermion interacting with a 17(1) gauge field. Gauge held 
in two dimensions is non-dynamical and can be explicitly integrated over. This generates 
a non-local four-fermionic interaction. Unlike the case of Thirring model, in this case we 
have a relevant interaction, and it affects the result for the entanglement entropy. 

It is well known that Schwinger model is equivalent to free massive scalar. We as¬ 
sume that single-interval entanglement entropy of Schwinger model maps to single-interval 
entanglement entropy of massive scalar held theory. It is not obvious that such an as¬ 
sumption is correct, because of non-local nature of duality between Schwinger model and 
massive scalar held [16]. Similar subtlety in relation between entanglement entropy in 
bosonic and fermionic pictures can appear even in CFT. It appears already in the case of 
multiple-interval entanglement entropy, which is not defined by simple expression analo¬ 
gous to (1.3), in CFTs of free boson and free fermion [24], Careful treatment shows that 
the results for multiple-interval entanglement entropy do indeed agree for the theories, 
related by bosonization [24], 

We match the expression for the entanglement entropy in Schwinger model with the 
result of [5] for massive scalar held. This gives an expression, which will be useful in 
calculation of entanglement entropy in the ’t Hooft model. We show that calculation of 
entanglement entropy in the’t Hooft model and the Schwinger model goes through similar 
steps. By explicit calculation we show that single-interval entanglement of un ction in the 
’t Hooft model (in the large -N limit) is given by a sum (with certain numerical coefficients) 
of entanglement c-functions of free scalars with masses in the spectrum of the’t Hooft 
model. 

The rest of this paper is organized as follows. In Section 2 we summarize the formulae 
for entanglement entropy, and outline the general calculation procedure which we will 
follow. We start with the CFT of free fermions, and then proceed to deforming the CFT 
by an arbitrary operator. We discuss in general how such an operator contributes to the 
single-interval entanglement entropy. In Section 3 we calculate entanglement entropy in 
the Thirring model, testing the formulae derived in Section 2. In Section 4 we discuss 
entanglement entropy in the Schwinger model, and derive expression, which will be useful 
for calculation in the’t Hooft model. In Section 5 we derive expression for entanglement 
entropy in the ’t Hooft model. We discuss our results in Section 6, where we provide an 
argument for why each meson, now matter how heavy it is, contributes to the entanglement 
entropy as a point-like particle. Appendix A is dedicated to derivation of gap equation 
governing dressed fermionic propagator in the ’t Hooft model. In Appendix B we review 
derivation of quark-antiquark scattering amplitude in the’t Hooft model. 
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2. Preliminaries 


In Subsection 2.1 we collect the defining formulae for entanglement entropy which we 
will need in this paper. In Subsection 2.2 we discuss free massless fermions on n-sheeted 
Riemann surface. In Subsection 2.3 we perturb the CFT of free fermions by operator O 
and discuss how it contributes to the entanglement entropy. 

2.1. Entanglement entropy 

Consider 2d system in a vacuum state. Suppose the subsystem A is a spatial interval 
of length L. Denote p to be density matrix obtained by tracing out all degrees of freedom 
outside of the A. One can compute the trace Tr p n for positive integer values of n. This 
quantity enters definition of the Renyi entropy 

S n = t —— log Tr p n . (2.1) 

1 — n 

Then, following the replica trick, one analytically continues the expression (2.1) to any 
real value of n, and calculates entanglement entropy using the formula 

S — — Tr p log p = — log Tr p n = lim S n . (2.2) 

on n-> 1 

n =1 

Evaluation of entanglement entropy therefore gets reformulated as a problem of finding 
the Renyi entropy (2.1). To calculate it, let us put the system on n-sheeted Riemann surface 
lZ n . Different sheets of lZ n are glued at the cut A. and boundary conditions across the 
cut are specified. Then the trace Tr p n can be determined as [1,2,3] 

Tr P n = |y (2.3) 

where Z n is partition function on lZ n . 

Due to (2.1), (2.3), the Renyi entropy can be determined as 

S n = —(log Z n - nlog Z 1 ) . (2.4) 

1 — n 

The entanglement entropy is therefore given by 

S = log Zi - ' ( 2 ' 5 ) 
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For the future purposes it is useful to express the partition function as 


^ n JJ2N Ai 


-Z(n ), 


( 2 . 6 ) 


where Ai = 0 and ^^\ n =i — 77 • Then entanglement entropy is 


0 N, T , ^ 1 dZ(n) 

S= T logL + logZ 1 -- — 


(2.7) 


n=1 


To get rid of the logarithm log Z\ in (2.7) one usually defines entanglement c-function 

( 2 . 8 ) 


_ d N cl 1 9Z(n) 

E dL 3 clL Z\ Bn 


n=1 


Suppose we start with conformal field theory with central charge equal to N. Then if 
we perturb it by some interaction, the entanglement entropy will be given by (2.7), where 


the y log L term is a CFT input [3], and Z(n) term originates from noil-trivial interaction. 


2.2. Free massless fermions 


Let us consider a free massless fermionic theory on n-sheeted Riemann surface 7 Z n . We 
have left-moving (holomorphic) fermion T = T i + bl , 2 and right-moving (anti-holomorphic) 
fermion T = Ti + with the free field action 

S=f d 2 w U*dV + W k dy) . (2.9) 

In our notation bar denotes chirality of the fermion, and star denotes complex conjugation. 
The integral in (2.9) is taken over n-sheeted Riemann surface lZ n with the branch points 
uq ; 2 , and the sheets are glued along the cut [wi,^]- Without loss of generality, we choose 
Wi — 0, W 2 — L, and take L to be real-valued. Instead of considering fermions (T, T) 
on lZ n , consider n fermions (T, , H/j), j = 1,... ,n, living on a complex 2 -plane [3]. The 
action on the complex plane is 


S = 



( 2 . 10 ) 


and we must insert twist operators T(0)T(L) into all correlation functions. Before describ¬ 
ing the action of the twist operators, let us do a conformal transformation of a complex 
plane, 

«=—, ( 2 . 11 ) 
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which maps twist operators to the points u — 0 and u — oo. These twist operators lead to 
the boundary conditions on fermions [4] 




—I —>• I — +1 

*j) \*j +1 




I:) - '- 1 * 


n+l 


Ti 


when one circles clockwise around u — 0, and 


\Ir 

T., 


->• 


T.-i 


Ti 


^(-1) 


n+l 




( 2 . 12 ) 


(2.13) 


when one circles around u = oo (he., when one circles around u — 0 counter-clockwise). 
The complex u-plane can be mapped onto a cylinder with world-sheet coordinates (r. a): 


u — 


e r+ia 


(2.14) 


The twist operators are inserted at r = ±oo, at infinite past and infinite future. The 
coordinate a G [0 , 27t] parametrizes circumference of the cylinder. 

Let us switch to diagonal basis of twist operators, by doing a discrete Fourier transform 


n — 1 
2 




e 27r ijs/n^ s ? 


Vi 


s 


n — 1 
2 

£ 

n — 1 

Z — 



e 2nijs/n^ 8 _ 


(2.15) 


In the basis ip s , ip s , s = — ..., the boundary conditions are 


(^ s \ k 2-Kis/n f ^ s \ 


(2.16) 


when one circles around u = 0 (goes around circumference of the cylinder), and 


f ^ s \ p -2nis/n ( ^ s \ 

wJ W,/ 


(2.17) 


when one circles around u — oo (goes around circumference of the cylinder in the opposite 
direction). The action for fermions (i/> s , ip s ) is 


S 


d 2 z 


n — 1 
2 

£ 


S = 


n — 1 
2 


ip*dili s + ip* s dil> s 


(2.18) 
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Fermions, satisfying the boundary conditions (2.16), (2.17), have the mode expansion 


^s( u ) = 

rez+ 1/2 

S(«) = 


h*») 



U^^ ’ 

, U r ~»+5 ' 
reZ+1/2 U 

r s | 1 5 

, U ™" r 2 

V’a(M) 

\ ^ ^00 

/ > _ r +iL + i 

U ' n' 2 

reZ+1/2 U 


(2.19) 


where the amplitudes of expansion satisfy anti-commutation relations [25,20] 

{V’(si) ri+Si /ni ^P(a 2 ) r 2 -s 2 /n} — ^s 1 ,s 2 ^r 1 ,-r 2 5 
ri — si/n’ ^( S2 ) r 2 +s 2 /n} ,s 2 ^ri, — r 2 

and the vacuum is defined as 4 

r+s/nl°) = 0 ’ 'lp(s) r+s/n\0) = 0 , r>0, 

(>) r+s/n I'0) = 0 , ^ (s) r+s/n |0) =0, r > 0 , 

for all s. Using equations (2.19), (2.20), (2.21) we compute the correlation functions 


( 2 . 20 ) 


( 2 . 21 ) 


(^* Sl (ui)^s 2 (u 2 )) = s. 

For the anti-holomorphic sector we obtain 

(VCiC(« 2 )) = s ai 


Sis 2 


S 2 


U2 

Ul 


Ul 

U 2 


U i — U 2 


( 2 . 22 ) 


Ul — U2 


(2.23) 


We can also arrive at the correlation functions (2.22), (2.23) by computing the limit 
of four-point functions of fermions and twist operators 


lim R 2As (a*(0)ip*(ui)ip s (u 2 )a a (R)) 

R —^oo 


1 


lim 


m — U2 -R-vxi y m — R ui 


U2 Ul — R 


U 1 — U2 \u 1 


U2 


(2.24) 


In (2.24) the A s = s 2 /(2n 2 ) is dimension of twist operator <J s (u) = e l ^^ s<<u \ We bosonized 
the fermions, ip s {u) = e l( ^ s ^ u \ where {^ s { u i)^s{ u 2 )) = — log (ui — M 2 ). Similarly one can 
derive correlation function for anti-holomorphic fermions (2.23). 


In our case s/n is never equal to 1/2, so we do not have to worry about zero modes: ^(s) 0 


never comes up. 



For the future purposes we also are going to need correlation function for fermions. 
Due to (2.15) and (2.22) we obtain 




1 e « u 


Uj ^ U^2 

2tt ij-^ _1_ 2'Kij 2 — 

e ri — e ri U% 


Similarly, for anti-holomorphic fermions we obtain 

rrij i 1 — X. 'Kijo 1 i 
,, . 1 e'Tiii" 2 e~~^~U2 n 2 

- (-i) J1+J2 --—- 2jr<J - 2 t 

/1 ^ - 777 n 


2irij 2 — 


e ri u i ~ c n U 2 


2.3. Interacting fermions 

Suppose CFT of free fermions is deformed by the interaction 


(2.25) 


(2.26) 


Sint = 9 d 2 zO(z,z ), 


(2.27) 


where O is a composite operator of T, T, and g is corresponding coupling constant. 
Suppose O in (2.27) is operator of dimensions (A, A). Then conformal transformation of 
O on complex plane, w = f(z), is given by 


8f(z)\ A fdf(z)\ A 


0(f(z),f(z)) = 0(z,z). 


(2.28) 


We can study the model (2.27) 011 Riemann surface 77.,, by considering n copies of the 
original model on complex 2 -plane, and specifying boundary conditions for going around 
singular points 2 = 0, L. Let us map the singular points to u = 0. 00 by transformation 
(2.11). The corresponding action is then 


Sint = g / d 2 u J ( u , u) Oj (u, u ), 


(2.29) 


and the boundary conditions are 


Oj(e 2ni u 1 e~ 2jri u) = O j+1 {u,u) , O n {e 2nl u, e~ 2ni u) = O x {u,u ). 


27 ri., ^ — 27ri—\ 


(2.30) 


In (2.29) we have introduced the J(u,u), which is a product of Jacobian of coordinate 
change (2.11) and factors coming from conformal transformation (2.28) of operators Oj, 


. L 2 - A - A (l-u) 2A (l-F) 2A 

J(M ' u) = -- 


(2.31) 
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Notice that boundary conditions (2.30) are just conditions for going between sheets of 
Riemann surface of the nth root function u 1 / n . Therefore the sum over j in (2.29) actually 
stands for summing over n sheets of u 1 ! n . This means that we can do the coordinate 
transformation u — v n , mapping n sheets of u 1 ! n onto complex plane with v coordinate. 
5 The way this transformation acts on O is described by equation analogous to (2.28), 

(nu n_1 ) A ( nv n - x ) A ^C>,(u n ,T n ) - 0{v,v). (2.32) 

j 

Now we can perform u = v n transformation in the action (2.29). Due to (2.32) we obtain 

Sint = gn 2 ~ A ~ A I d 2 v J(u n ,T n ) (V _A u 1_A ) 0(v,v). (2.33) 

Up to this moment the discussion has been general and did not rely on the fact that 
Oj operators are composed of T ;/ , fermions. After we do the change of variables u — v n 
we obtain operator 0(v,v) in (2.33), composed of single-valued on a plane fermions T(u), 
T(tJ). The free-held correlation functions of these fermions are 

<**K)*(t>2)> = ——, <T(?7 1 )*(v j )) = ^2_. (2,34) 

Vi —V2 Vi — V 2 

This conclusion can also be reached in a more explicit way. First, one writes down 
perturbative expansion of the partition function 

Z = J [DT] e ~ Sint , (2.35) 

with Si n t given by (2.29). Each term in the perturbative expansion can be evaluated, by 
using correlation functions of fermions (2.25), (2.26). One notices that it is possible to do 
the transformation u — v n for each given order of the perturbative expansion. Indeed, for 
example in holomorphic sector every «m" is multiplied by phase factor e 27r *im,/n^ aric j j m 
is summed over from 1 to n, see (2.25). Then all the terms can be assembled back into 
partition function of fermions (2.34), deformed by interaction (2.33). 

5 Mapping from 7 Z n to complex plane is frequently used in calculations of entanglement entropy, 
see, e.g., [3,26]. 
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3. Thirring model 


As a warm-up, in this section we consider Thirring model, described by the action 

S T — J + A T W*WWV) , (3.1) 

where At is a dimensionless coupling constant. A well-known fact is that one can imme¬ 
diately bosonize the fermions, W(z) = e l v( z \ T(;7) = and notice that the action 

(3.1) describes free massless scalar field <p. Therefore the entanglement entropy of a sin¬ 
gle interval in Thirring model is simply given by (1.3) with c = 1. We prefer to postpone 
bosonization and instead treat the four-fermionic term in the action (3.1) as a perturbation 
O. Then we can test the procedure, described in Subsection 2.3. 

Let us put the model (3.1) on n-sheeted Riemann surface 7 Z n . This can be accom¬ 
plished by dealing with n copies of fermions, (Tj, Tj), j — 1 ,..., n on a complex z-plane, 
with the action 


/ n 

d 2 z^ (Vjd'&j + A / , 

j =i V 


(3.2) 


and twist operators inserted at z — 0, L. The interaction in (3.2) is of the form (2.27) 
with (A, A) = (1,1). In the u coordinates (2.11), the fermions (Tj, Tj) satisfy boundary 
conditions (2.12), (2.13). 

Due to (2.31) we obtain J(u,u ) = 1. The partition function of Thirring model is 
therefore given by 


Zi = 


L 2A - 


where 


[DV]e 


A n - 


-A T f d*u'£,"** j *j*** j 


n 2 — 1 
12 n 


(3.3) 


(3.4) 


We can apply conformal perturbation theory for computation of (3.3). In the previous 
section we computed CFT correlation functions (2.25), (2.26) for the (Tj, T ? ) fermions. 
We can either use these correlation functions for explicit calculation, or use general result 
derived in Subsection 2.3. Due to (2.33) we obtain simply 


7 1 
L 2A - 


^ —A t J d 2 vdipd(p\ 


/ ’ 


(3.5) 


where we have bosonized the fermions, and 


(^(ui)^(u 2 )) = -log(vi - v 2 ), (<p(ui)<p(u 2 )) = -log(m - v 2 ) ■ (3.6) 

From (3.5) we conclude that entanglement entropy of Thirring model is given by CFT 
formula (1.3), with c = 1, as expected. 
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4. Schwinger model 

In this section we are going to derive formula for entanglement c-function in the 
Schwinger model (2 d QED with massless fermion). Derivation splits into three parts. 
First, by using the procedure described in Subsection 2.3, we derive an expression for Z(n ), 
defined in (2.6), which contains all the lion-trivial contribution from the gauge interaction. 
The resulting formula contains fermionic four-point function in Schwinger model, which 
we compute in Subsection 4.2. We assemble all the results in Subsection 4.3 to write down 
expression for entanglement c-function. 

The purpose of this section is twofold. First, it can be viewed as a warm-up for 
’t Hooft model: in this section we study entanglement entropy in 17(1) gauge theory with 
massless fermion, and in the next section we generalize it to U(N ) gauge theory with 
massless fundamental fermion. Second, Schwinger model is equivalent to a free massive 
2d scalar, for which entanglement entropy is known (numerically) [5]. It turns out that 
calculation of entanglement entropy in Schwinger model is very similar to calculation in 
7 Hooft model, as we will show in the next Section. Therefore the result for 7 Hooft model 
can be expressed using known result for Schwinger model. 

4-1- Derivation of partition function Z m (n) 

The Schwinger model in the 7 = 0 gauge, with the gauge coupling m, is described by 
the action 

S = j d 2 z ^-^(<97) 2 + m7T*T + . (4.1) 

The gauge field is non-dynamical in two dimensions, so we can integrate over it. This 
amounts to substitution of 

A(zi,zi) = m J d 2 z 2 G{z 1 - z 2 )^*{z 2 )^{z 2 ), (4.2) 

into the action (4.1), where 

d 2 G{z\ — z 2 ) — 5{z\ — z 2 ) =>• G(zi — z 2 ) = ——-— • (4-3) 

2-7T Zi — z 2 

Therefore we can consider Schwinger model as a CFT of free fermions perturbed by the 
interaction of holomorphic fermions 

Ss = m 2 j d 2 z 1)2 G(zi - z 2 )^(z 1 )^{z l )^{z 2 )^{z 2 ). (4.4) 
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It is well known that the Schwinger model is equivalent to a free scalar with mass m. This 
can easily be seen by bosonizaion the fermions. Then (4.4) becomes just a mass term for 
a free scalar. As in the case of Thirring model, let us postpone bosonization for now. 

Let us consider Schwinger model on n-sheeted Riemann surface 7 Z n , with the sheets 
glued along the cut of the length L. As in the case of Thirring model, this can be ac¬ 
complished by considering n copies of fermions (47j, ) on a complex plane, with the 
interaction given by 

/ n 

d 2 z lj2 G{z i - z 2 ) 22 • (4-5) 

j-1 

By doing change of variables (2.11) we obtain 

/ n 

d 2 u lt 2 F{u r, u 2 ) 22 , (4.6) 

3 = 1 

where we have denoted 

F(ui,u 2 ) = - -- t=G(ui — u 2 ). (4.7) 

|1 — Ui\ z \\ — u 2 \ z 

The action (4.6) can be written in a form, similar to (2.29), 

/ n 

d 2 u lj2 F(m, u 2 ) y j Oj{ui)Oj{u 2 ), (4.8) 

5 =i 

where operators Oj(u) = 'L*(n)'I' ) (n) of dimension (A, A) = (1,0) satisfy boundary condi¬ 
tions (2.30). Therefore for the model (4.8) to be well-defined, the coordinates u± }2 should 
go around singular points u — 0, oo at the same time. This means that the sum over j 
allows to make change of variables Ui. 2 = Similarly to (2.32) we obtain 

ArV 1 H W) = OMOM , (4.9) 

3 

where 0{y) = 4/*(n)\l/(n), and correlation function for T(n) is given by (2.34). The action 
therefore becomes 

S s = n 2 m 2 L 2 J d 2 v ia H(v u v 2 ; n)T*(ni)T(m)T*(n 2 )T(n 2 ), (4.10) 
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where we denoted 


—n — 1 —n — 1 


H(vi,v 2 ;n) = 


II 


,n|2 


1 


,n\2 


(4.11) 


J 1 I I ^ — u 2 I — u 2 

Again, one can also arrive at the same result (4.10) by using correlation functions (2.25), 
(2.26) in explicit perturbative calculation of the partition function. 

The partition function of Schwinger model on 7 Z n is given by 


Zn ~ L 2A n 


[.DT] e 


—Ss 


(4.12) 


In (4.12) we have path integral in CFT of free fermion \I/(u) with correlation function 
(2.34), perturbed by interaction (4.10). Following (2.6), let us denote 


Z m {n) = / [T>T] e 


— Ss 


(4.13) 


Therefore Z{ 1) = Zf is just a partition function of free massive 2 d scalar on a complex 
plane. For the purpose of calculation of entanglement entropy (2.7), we are interested in 

= —m 2 L 2 f d 2 v 1 , 2 H 0 (v 1: v 2 ) Zf ^*(vi)tf (ui)tf*(v 2 )tf (u 2 ) e - Ss| "=^ , 

n= 1 

(4,14) 

where we have denoted 


d Z m (n) 
dn 


H 0 (vi,v 2 ) = 2 i7(m,n 2 ; 1) + 


dH(v!,v 2 ;n) 


dn 


(4.15) 


71=1 


Now we need to compute fermionic four-point function, which enters (4.14). 


4-2. Calculation of quark four-point correlation function 

To compute correlation function in (4.14), let us bosonize the fermions, 

— rn 2 L 2 J cl 2 vi j2 Hq(vi 1 v 2 ) (di<&{y\,vi)d 2 <&(v 2 ,v 2 ) e~ s<s> ) , (4.16) 

where due to (4.10), (4.11) we have 

S$ = —m 2 L 2 [ d 2 v x , 2 T7-iwr:- jozr -9i$(m, vi)d 2 ®(v 2 , v 2 ) . (4.17) 

J |1 - Vi\ |1 - v 2 \ vi ~ v 2 


1 dZ m (n) 
Zf dn 
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We can easily get rid of the extra factors in the integral (4.17) by doing transformation, 
inverse to (2.11): 


2 /1,2 ~ 1 
2 / 1,2 


= W,2 , 


(4.18) 


where y\^ are new dimensionless coordinates on complex plane. This gives 


2 t 2 


S$ = —mL 


,2 y 1 2/2 0 T 2 r 2 

d 2 / 1,2 =-— = m L 

Vi ~ 2/2 


d 2 y $(2/, 2/) 2 , 


(4.19) 


where we have also used (4.3). Applying the same transformation in (4.16) gives 


1 d Z{n) 
Zf dn 


2 r 2 


m L 


d 2 yi,2 Hi(yi, y2)d yi d y2 ($(2/1,2/1) $(2/2,2/2) e S *> 1 


(4.20) 


where we have denoted 


= Ho_T. (4.21) 

V 2/1 2/2 / 2/12/2 

Let us perform Wick rotation and go to light-cone coordinates, (y. y) — (y + . y~), 
( p,p ) = (p+,p-). Then p ■ y = pyy + + p~y~ = py + py■ We are working in A- = 0 
gauge. Therefore the residual gauge transformations are parametrized by a = a(y + ). In 
the light-cone (y + ,y~) coordinates it becomes manifest why (4.14) is a gauge-invariant 
object, since fermions as not separated in the y + direction, yf — y 2 — y + ■ 

Due to (4.19) we have in (4.20) correlation function of free massive scalar, with di¬ 
mensionless mass m 2 L 2 , 

(•Hy + ,yT)t(y+,yi)e- s °)= j r-» 3 -) _L_ . ( 4 , 22 ) 

Therefore, from (4.20) we arrive at 

= rn 2 L 2 I dy^ 2 dy + H 1 (y 1 ,y 2 ) J d 2 qe iq ~ [y i ~ y a } ^~ m 2 L 2 • ( 4 - 23 ) 

The integral (4.23) will also appear in the calculation of entanglement entropy in’t Hooft 
model. This integral requires regularization, and as written, should be viewed formally. In 
the next Subsection we derive a useful expression for this integral. 


1 dZ m [n ) 
Zf dn 
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4-3. Entanglement c-function 

Due to (2.7) to find entanglement entropy we need to evaluate (4.23). This is a hard 
exercise and, as mentioned above, requires regularization. Fortunately we know the answer 
for entanglement entropy S m (mL ) of free scalar field with mass m [5,6,18]. Therefore we 
can use the known value of S m (mL) to find (2.7) 


1 d Z rn (n) 


Z f 


dn 


= \ log L + log Zi - S m {mL). 


(4.24) 


n=1 


Differentiating w.r.t. L we obtain 


dL 


1 dZ m (n) 


Zf dn 

n=1 - 


= i - c m (mL ), 


(4.25) 


where c m (mL) is entanglement c-function of free scalar with mass m. Equations (4.23), 
(4.25) give rise to 


dL 


m 2 L 2 / dy 12 dy + H 1 (y 1 ,y 2 ) I d 2 qe 


' 2 np iq - ( ' y i ~ y 2 ). 


Q- 


q 2 — m 2 L 2 


= i - c m (mL) 


(4.26) 

Equation (4.26) is the main result of this section, and it will play a key role in derivation of 
the entanglement entropy in’t Hooft model. When m — 0 we have Co = 1/3, and therefore 
regularized value of the expression in the l.li.s. of (4.26) is equal to zero. 


5. ’t Hooft model 

In this section we derive expression for entanglement of un ction in the’t Hooft model. 
As in the case of Schwinger model, derivation splits into three parts: calculation of Z(n), 
defined in (2.6); calculation of fermionic four-point function in the ’t Hooft model; and 
assembling all the results together, to get the expression for the entanglement c-function. 
Calculation of Z(n) is similar to analogous calculation, performed for the Schwinger model. 

5.1. Derivation of partition function Z[n) 

The’t Hooft model in the A" = 0 gauge is described by the action [11] 

S = f d 2 z (~{dA a ) 2 + + g1‘A< 0 . (5.1) 
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Here T a , a — 1,..., N 2 are generators of the u(N) algebra, and g is coupling constant, with 

- Oi 

dimension of mass. Fermions \E ,a , T , a = 1,..., N live in fundamental representation of 
the U(N'). Similarly as in the case of Schwinger model, we can integrate over the gauge 
field, which amounts to substitution of 

A a ( Zl )=gJ d 2 z 2 G( Zl - z 2 )V a *(z 2 )T^(z 2 ) (5.2) 

into the action (5.1). This gives rise to free fermionic theory, deformed by the interaction 

S H =g 2 J d 2 z h2 G( Zl - z 2 )^ a *( Zl )^( Zl )^( Z2 )\t a (z 2 ) . (5.3) 


In (5.3) we have used 


t; 


a 


-i a 


aiO!2 Q!3Q!4 


=u 


Oii Q!4 V(X2OL3 


(5.4) 


The Green’s function for d 2 operator is given by (4.3). 

Let us consider’t Hooft model on n-sheeted Riemann surface 7 Z n with branch points 
at z = 0, L. Let us make coordinate change (2.11), and consider n copies of fermions 
(T", 47 ? ) on a u-plane, satisfying boundary conditions (2.12), (2.13). These fermions are 
described by CFT, deformed by the interaction 


S H = g 2 L 2 J d 2 u t ^ F(u u u 2 ) £ , (5.5) 


where F(ui,u 2 ) is defined in (4.7). 

Let us proceed as in the case of Schwinger model. The action (5.5) can be written as 


S H — g 2 L 2 J d 2 u 1 , 2 F(u 1 ,u 2 )J2of(.'“i)°j a M, (5.6) 

j 

where (u) = 47“*('Ui)47^('Ui). For this action to be well-defined the u\ )2 should go 
around singular points Ui j2 = 0, 00 at the same time. This means that we are allowed 
to make a change variables Ui j2 — Vy 2 • mapping n sheets of u 1 / n function onto complex 
u-plane. The transformation law of operators O a & is given by 


n 2 v^ x vl~ x ^ Of(v?)0? a (y%) = O aP { Vl )Od a (v 2 ) , ( 5 . 7 ) 

j 

where O a/3 (v) — 47“*(u)\I7^(u). Correlation function of holomorphic fermions on w-plane is 

(^“> 1 )^ 2 )) = 5 ap ---. ( 5 . 8 ) 

vi - v 2 
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In the v coordinates the action (5.6) becomes 


S H = n 2 g 2 L 2 f d 2 v 1}2 H{ Vll v 2 ; n)^ a *( Vl )^(v^*(v 2 )^ a (v 2 ), (5.9) 

where H(vi,v 2 ;n) is defined in (4.11). 

The partition function of’t Hooft model on 7 Z n is therefore given by 

Z " = /P*] ^ S " ■ (5-10) 

In (5.10) we have path integral in CFT of free fermions with correlation functions (5.8), 
perturbed by interaction (5.9). As in the case of Schwinger model, and following (2.6), let 
us denote 

Z{n) = j[D^]e~ SH (5.11) 

in (5.10). To compute entanglement entropy (2.7), we have to find 

= -g 2 L 2 j d 2 v h2 H 0 (v u v 2 )Z? , 

(5.12) 

where Hq(vi,v 2 ) is defined in (4.15), and Z^ 1 is partition function of’t Hooft model on a 
plane. In (5.12) we have SA| n= i, which is just four-fermionic interaction in’t Hooft model 
on complex plane. Indeed, doing change of coordinates (4.18), we obtain 

= -g 2 L 2 f , (5.13) 

n= 1 

where H\(y\. y 2 ) is defined in (4.21). Subscript H stands for correlation function in the 
theory with interaction 

Sh = g 2 L 2 [ d 2 Vl , 2 W a *(y 1 )^(y 1 )^*(y 2 )^“(y 2 ). (5.14) 

J Vi - y 2 

We can bring the gauge field back, restoring the original’t Hooft model framework, with 
fermions interacting locally with the gauge field: 

sh=j i 2 y (y-\{0A a ? + gA° , (5.15) 

where 

g = gL (5.16) 

is dimensionless gauge coupling. 


1 d Z(n) 
Z[ r dn 


d Z(n) 
dn 
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5.2. Calculation of quark four-point correlation function 


In (5.13) we have four-quark amplitude in’t Hooft model on a plane, 

K.(y i.rs) = (y a *(yi)'H ,3 (.yi)'l p *(y2)y a (.y2)) H ■ (5.17) 

Due to translational symmetry, £(yi,y 2 ) = £( 2 /i — 1 / 2 )- This is gauge-invariant quantity, 
which has been demonstrated explicitly in [27]. The situation is analogous to Schwinger 
model, discussed in the previous Section. Remember that we are working in the gauge 
A a = 0, and therefore the residual gauge transformations are parametrized by a a such 
that 

0 = 5A a = da a + f abc A b a c = da a , (5.18) 

that is, a a — et a (z). Therefore the object T q *( 2 /i)T“( 1 / 2 ) is gauge-invariant simply because 
both fermions are not separated in the y direction. Therefore the four-fermionic amplitude 
in (5.13) is gauge-invariant. In fact one can consider a simpler gauge-covariant “blob” [11] 

S^ a ^ r ( yi ,y 2 ) = (r|T Ql *( yi )T“ 2 (y 2 )|0>, (5.19) 

giving an amplitude of creation of meson |r) by quark-anti-quark pair T a *(i/i)\I/ a (i/ 2 ). 
The Fourier transform of amplitude (5.17) is given by 

JC(yi-y 2 )= I rfV,2,3,4e i(?/1 ’ (Pl - p2)+J/2 ' (P3 - p4)) (* a * ( Pl )& (p 2 )§^ (p 3 )V a (p 4 )} r . 

(5.20) 

Define vertex T(p,p';q ) [12] for quark-anti-quark scattering as 

^ ai *( Pl )^ a2 ( P 2)^ a3 \psW 4 (p4)) H =S(p 1 -p 2 +P3-P4) P(Pi)P(P2)P(P3)P(P4) 
x T( P1 ,P 2 ; P1 - P4 ) ~ S aia2 S asa4 T( Pl ,p 4 ; Pl - p 2 )} , 

(5.21) 

where P(p) is the full quark propagator [11,12], Let us use (5.21) in (5.20), obtaining 

£( 2/1 - y 2 ) = £ 1 ( 2/1 - y 2 ) + £ 2 ( 2/1 - 2 / 2 ), (5.22) 


where 

£ 1 ( 2/1 - 2 / 2 ) =N 2 J d 2 p h 2 A e i ^-^<y^T(p 1 :P 2-,p 1 -p 4 ) 
X P(j)l)P(p2)P(-Pl +P2 +P4)P(P4) , 


(5.23) 
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£2(2/1 - 2/2 ) = N 


d 2 p i, 2 , 4 e <(pl_P2) ' (!/1 " w) 


T{Pi,PAPi -P2) 


(5.24) 


X P(p 1 )P(p 2 )P(-pi +p 2 +P4)P(P±) , 

The iV 2 in /Ci appeared after we summed over op = 04 = a and a 2 — «3 = (3, and AC in 
AC 2 appeared after we summed over a — p. 

Let us denote integrated momenta in (5.23) as pi — p, p 2 — p'. P 4 — p — q, and define 
Ki(q)= I d 2 pd 2 p e i{p - p,) < yi - y2) T(p 1 p'-q)P{p - q)P(p)P{p - q)P(p). (5.25) 


Therefore 

£ 1 ( 2/1 - y 2 ) = N 2 J d 2 qK 1 (q ). (5.26) 

Let us denote integrated momenta in (5.24) and pi = p, p 2 = p — q, P 4 — p\ and 
dehne 

K 2 (q) = J d 2 pd 2 p' T(p,p';q)P(p - q)P(p)P(p' - q)P(p '). (5-27) 

Therefore 

£2(2/1-2/2) = N J d 2 qe iq ^~ y2) K 2 (q). ( 5 . 28 ) 


The AC 2 is sub-leading in the large-iV limit, and we will not consider it. In what follows we 
will only be interested in ACi. 

Let us perform Wick rotation and switch to light-cone coordinates, similarly to what 
we have done in Subsection 4.2, for Schwinger model. Denote g to be gauge coupling 
constant. The ’t Hooft coupling constant is A = r/A r2 . The extra hat over the couplings 
is introduced because at the end we will substitute g — gL , where g is the original gauge 
coupling constant, and L is the length of interval, see (5.16) and preceding discussion. 

Let us introduce the IR cutoff A/ r. Then full quark propagator, in the large-AC limit, 
is given by [ 11 ] 6 


Pip) 


2 _ Asgn(p^) 

1 ^ 7T A TR 


M 2 —ie ’ 
V- 


(5.29) 


where 

M 2 = m 2 — A/ 7 r, 


(5.30) 


6 Equation for full quark propagator P(p) is derived in Appendix A. Quark propagator is two 
by two matrix. In A- — 0 gauge we are only considering entry of this matrix which corresponds 
to two-point function of left-moving quarks, since right-moving quarks are non-dynamical. 


20 



and fh q — m q L is dimensionless bare quark mass. At the very end of calculation of gauge- 
invariant observables one sends A//,> —> 0. 

It is known that T(p,p';q ), in the large -N limit, satisfies equation [12] (see Appendix 
B, where derivation of this equation is reviewed) 


T(p,p';q) = 


ig 


(p- — p'_) 2 7 r 2 


iX f „ 1 

H-o / dk— 


( k _ — p-) 2 


dk + P(k)P(k-q)T(k,p';q). (5.31) 


Denote 


P- 


/ P- 
x = — 


x — z — , x' = — , (5.32) 

q- q- 

which stand for the fraction of total momentum g_, carried by one of in-coming quarks 
P-, and one of the out-going quarks p'_. The solution to equation (5.31) is [12] 


T(p~,p'-', q) = 


A 


A 2 x ^ 1 

n 2 _ 


N(p_ — p'_) 2 7rNq_ q z — frif, + ie j 0 


, f 1 €(y')My) 

Jo (</' -1') 2 to - s) 2 ' 

(5.33) 

Integrals over y and y' in (5.33) are assumed to be taken with the cutoff: holes of size 
A m/q~ are drilled around y — x and y' — x'. In (5.33) we denoted (j> r (y) to be set of 
eigenfunctions of’t Hooft equation 


7rM ‘ 2 a r \ td f 1 a ' Mv') M 

x^V) My) - p Jo dy ¥^W = T Myh 


(5.34) 


describing meson of mass m r — m r L. Here P stands for principal value. 

Eigenfunctions <f> r (x) of’t Hooft equation (5.34) are non-vanishing when x G [0, 1] 
and satisfy completeness and orthonormality conditions [11,12] 


cj) r (x)(j) r (x') = 5(x — x '), 


dx (j) r (x)(f)*,(x) = S r 


(5.35) 


For large r (since N does not appear in (5.34), large r is not 0(N P )) we have [11,12] 

4> r (x) ~ v^sin (itrx ), m 2 ~7r 2 Ar. (5.36) 


Let us re-write (5.33) as [21] 

A 4A 2 


T(p-,p'_;q) = 


+ 




N(p- — p'_) 2 7lNX 2 IR g2 _ f n r 2 _J_ 


Xr{x,qJ)Xr{x',qJ), (5.37) 
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where we have defined 




(5.38) 


The reason for such reformulation is that one can see that presence of A / R in the r.li.s. of 
(5.38) makes the Xr finite in A ir —> 0 limit [21], Indeed, we can split the integral over 
y into principal value part, which is finite when Xjr —>> 0, and infinite part, coming from 
divergence near y — x: 


We conclude that 


L dv W=W = P L + 


Xt(x, Q-) = 4>r{x) + 0(X IR ) . 


(5.39) 


(5.40) 


Now let us discuss /Ci(yi — 1 / 2 ), see (5.25), (5.26). We have 


£i(g) = / dp- dp'_ e^-p'-Hv, -v 2 ) T(p_,p'_; q) h I 2 . 


(5.41) 


we have introduced integrals I\ . 2 , 


Ii — J dp + e ip+ '( y i P(p — q)P(p ), 

/.. = / dp' e-v'^vi-vi) /•(,/ _ q )P(p'). 


(5.42) 


As discussed above, here we should set y{~ — y 2 ■ In A i R —* 0 limit 


h = h = 


m 2 \iR 


(5.43) 


Using definition (5.32) and equations (5.37), (5.40), (5.41), (5.43), and neglecting overall 
factors of 7 r, etc., we obtain 


1 /*1 _ _ q 2 

KAq) = — / dxdx'e i{ - x ~ x ' )q -^ ~ y i } 

N J 0 “ q‘ ~ 


(j) r {x)4) k r {x' 


(5.44) 


Now we need to plug (5.44) into (5.26) and integrate over q, 


'' £ 1 ( 2/1 - 2 / 2 ) = A f dxdx' j d 2 q e^ x x y 2 ) ^ 7 ,<f> r (x)<f)*(x') (5.45) 


g* — rrq 
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We can change variables of integration (assuming some smooth regularization when x = x', 
which we discuss below), g_ —» q~/(x — x '), q + —* (a; — a/)^. Then (5.45) becomes 


/Ci( 2/1 ~y 2 ) = A / dxdx' 


(x — x 


- 2 Jd?qe*-K-MYi 


q 2 — m 2 


<f>r(x) 4>r( x ') ■ 


(5.46) 

Divergence in the integral in (5.46) when x approaches x' originates from divergence in the 
IR limit, when momentum —p'_, flowing between ij\ and 1 / 2 , goes to zero. Therefore we 
can assume that such integral is regularized as in ’t Hooft equation, by taking principal 
value part: 


dx dx 


(j) r (x) 


' 7 - -K(j)r(x)(j)*(x , ) = - dx (j)*(x') P / dx -—-—. 

[x-x'Y J o J o ( x-x'Y 


(5.47) 


Using’t Hooft equation (5.34), normalization condition (5.35) of’t Hooft wave-functions 
4> r (x), and definition (5.30), we can re-write (5.47) as 


mf / , 

Cl r — “^—h 1 

A V 


A Jo 


dx \(j) r (x)Y 


x(l — x) 


(5.48) 


The integral in (5.48) is positive-valued, and for large r we can use (5.36) and substitute 


dx \J) r {x)\ 2 — -- ~ log(27rr). 

x(l — x) 


(5.49) 


In this subsection we have been considering the situation of generally non-vanishing 
bare quark mass rn q . However, our computation of entanglement entropy is performed in a 
set-up when we perturb CFT by four-fermionic interaction (5.9). If we had also perturbed 
CFT by non-vanishing bare mass for quarks, then an extra term, besides (5.9), would 
have been generated, and it would have also modified the total contribution (5.12) to the 
entanglement entropy. In particular, it is known that entanglement entropy of free massive 
fermion is given by a non-trivial expression [4,6]. Therefore let us refrain to vanishing bare 
quark mass. Then from (5.48) we obtain that a r > 0 for all r, 


a r = + / dx \4> r {x)\‘ 


A Jo 


x(l — x) 


(5.50) 


5.3. Entanglement c-function 


We can use the result (5.46) in the expression (5.13), 


1 dZ(n) 
ZY dn 


a a r J d yi,2 d v +H Yyi,y2) J d 


2 qe ^-^-y^—^— 

q z — rn z 


(5.51) 
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To find entanglement c-function in the ’t Hooft model we need to substitute (5.51) into 
(2.8), and compute the integral. This is lion-trivial, and requires regularization. In fact, 
the integral which appears in (5.51) is exactly the same as the integral which appears in 
analogous expression for the Schwinger model (4.23). 

To avoid the difficulty of calculation of integral in (5.51), we use (4.26) to obtain 
entanglement c-function in the’t Hooft model, C//(A), in terms of entanglement c-function 
of free massive scalars c(m r ) 


Cff(A) = 



(5.52) 


where m r = m r L. We know that Cij(X) should be a monotonically decreasing function 
[18], which should fix unambiguously overall sign of the sum in (5.52) to be plus (at the 
end of previous Subsection it was shown that a r > 0, see (5.50)). Due to (5.50), (5.36), 
(5.49), we have a r X/fh 2 = 1 + 0(r~ l log r) for large r. 

The expression (5.52) cannot be a final answer, and requires regularization. This can 
be consequence of the way we handled integrals in (4.26) and (5.51). Besides, as written 
now, (5.52) has a contribution TV/3 from twist operators, which diverges in the large-TV 
limit, and the sum over mesonic states, which is finite in the sense of large-TV counting. 
To fix this problem, let us take non-trivial L-dependent part of (5.52), and regularize the 
sum by requiring correct infra-red behavior. Namely, we have ch\l=oo — 1/3, which is the 
entanglement c-function for one free massless meson (mi = 0) in the spectrum. This way 
we arrive at 

1 °° A 

c//(AL 2 ) = - + V] a r ^ c(m r ) • (5.53) 

3 mi 

r=2 

Recall that [5,6] entanglement c-function of free massive scalar falls off exponentially at 
large distances, c(m) ~ e~ m . The masses of higher mesons are given by (5.36). Therefore 
for interval of length L the number of mesons which contribute to the sum (5.53) is 


O 



(5.54) 


As L is decreased, and we go to the ultra-violet, more mesons start to contribute, and 
when XL 2 = 0(1/N) we get Ch — TV due to (5.54), with 0(N ) light mesons contributing. 
In this regime of parametrically small ’t Hooft coupling constant we can neglect gauge 
interaction of quarks. Therefore we need to switch to picture of free quarks, which gives 
ch — TV/3. It would be interesting to interpolate more smoothly between the regimes of 
free quarks and mesons. 
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6. Discussion 


In this paper we used the replica trick approach to calculate entanglement c-function 
in three models of a two-dimensional quantum field theory. We have considered Thirring 
model, Schwinger model, and ’t Hooft model. Thirring model is the simplest model with 
four-fermionic interaction, and after bosonization it becomes re-formulated as a theory 
of free massless scalar. Postponing bosonization, we have demonstrated explicitly that 
single-interval entanglement entropy in the Thirring model is given by the CFT expression 
(1.3) with central charge c = 1. 

We have shown that replica trick calculation of entanglement entropy in the Schwinger 
and’t Hooft models follows a similar path. In fact, one can derive entanglement entropy in 
the ’t Hooft model, provided it is known what is the entanglement entropy in the Schwinger 
model. The latter is indeed known (numerically in general and analytically in asymptotic 
regions), because the Schwinger model is equivalent to free massive scalar field. 

It is essential for our calculation that ’t Hooft model is solvable in the large -N limit. 
More precisely, it is important that one can compute quark-anti-quark scattering ampli¬ 
tude, and demonstrate that the only intermediate states in this process are mesons. Since 
each meson is free in the large-fV limit, we can use the known result for entanglement 
entropy of free massive 2d scalar. We have demonstrated by explicit calculation that en¬ 
tanglement c-function in the ’t Hooft model (in the large-IV limit) is equal to a sum of 
entanglement c-functions of all the mesons in the spectrum (with certain numerical coef¬ 
ficients). It is interesting to observe that the fact that each meson has a finite size, which 
moreover grows as the meson gets heavier, is not manifested in the entanglement entropy. 

A possible explanation of the fact that each meson contributes to the entanglement 
entropy as a point-like particle is the following. As it was mentioned in Introduction, 
defining entanglement entropy in systems with gauge fields is subtle. The problem appears 
with finite-size gauge-invariant operators (like mesons in the ’t Hooft model) which are 
partly located in A and partly in A. Such operators are in conflict with the very definition 
of entanglement entropy, which starts with the assumption of splitting of the Hilbert space 
into a tensor product 

n A xn-j. (6,i) 

Therefore the question is what exactly the replica trick calculation, which we performed 
in this paper, has given us. 


25 



One answer is that in the entanglement entropy calculation the gauge-invariant op¬ 
erators, which are located partly in A and partly in A. are cut into two parts, by the 
entangling surface. The free ends of those parts are glued to the extra added gauge non¬ 
invariant (charged) edge inodes, located on the entangling surface [16,28] (see also [29] 
for the calculation of the contribution of the edge modes to the entanglement entropy). 
Therefore each of the resulting parts is completed to a gauge-invariant object. One is 
located entirely in A , the other is located entirely in A. This procedure is incorporated in 
the replica trick calculation [16]. This means that the whole Hilbert space is enhanced to 
include such charged boundary degrees of freedom, located on the entangling surface. The 
consequence of such an enhancement is that there are no more gauge-invariant operators, 
located partly in A and partly in A. This makes definition of the entanglement entropy 
consistent. At the same time it explains why finite size of mesons is not manifested in the 
entanglement c-function, no matter how long the meson is. 7 

If bare quark mass is zero, then the lowest meson is massless. Therefore in the deep 
infra-red we have CFT of one free massless scalar, and entanglement c-function is equal 
to 1/3. In the deep ultra-violet regime quarks are free, and entanglement c-function is 
equal to N/ 3, where N is the number of colors. One can compare this crossover from the 
mesonic phase to free quark phase with the phase transition, observed in [31] (see also 
[32]), where single-interval (slab) entanglement entropy in confining gauge theory has been 
calculated holographically. In [31] it was shown that there is a critical thickness of the 
slab L c so that when L < L c entanglement entropy behaves as N 2 , and when L > L c 
entanglement entropy behaves as A r0 . It would be interesting to understand better how 
these two observations are related. 

It would be interesting to address the question, raised in the introduction, and calcu¬ 
late entanglement entropy in the quark language, perturbatively in the coupling constant. 
Such calculation will shed light on the crossover region, XL 2 — 1/N. 

1 In a pure two-dimensional Yang-Mills theory there is an additive contribution to the entan¬ 
glement entropy, which can be interpreted by counting the gauge non-invariant states, located on 
the entangling surface [16]. See also [30], where entanglement entropy in the 0(N ) cr-model has 
been calculated, and interpretation in terms of counting the states of the UV degrees of freedom 
has been provided. It might be that there are such parton-counting terms, entering entanglement 
entropy in the ’t Hooft model, and these terms are L- independent, and hence invisible in the 
calculation of the entanglement c-function. We thank A. Wall for discussion of this point. 
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Appendix A. Gap equation for fermionic propagator 

In this Appendix we derive gap equation on full quark propagator by calculating path 
integral in the large-A limit. This is different from the conventional derivation, which 
involves equation satisfied by planar ladder diagrams [11], Similar derivation for Gross- 
Neveu model can be found in [33,34], 

For generality, let us give fermions the bare mass m. The total action is then 

S = S 0 + S int , (A.l) 


where 




cl 2 z (2)(7 t ‘3„ + m)$“(2) 


(A.2) 


and (let us absorb ’t Hooft coupling A into rescaling of coordinates Z \, 2 and the mass m) 


S in t = ^ J d 2 z ia G( Zl - z 2 )^ a *(z 1 )^ a (z 2 )^(z 2 )^(z 1 ) . (A.3) 


We have denoted Dirac fermion as T, and its charge conjugate as 4/t, 



(A.4) 


and the Dirac operator is 



(A.5) 


Let us introduce auxiliary gauge-singlet field P(z\,z 2 ), which for the purpose of finding 
quark two-point function, due to translational symmetry, is taken to be P(z\ — z 2 ). We 
can re-write the action (A.3) as 


Sint = I d 2 z 1)2 G{ Zl - z 2 ) (NP(z\ - z 2 )P(z 2 - zi ) + P{z x - z 2 )^ a *( Zl )^ a (z 2 )) (A.6) 
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Then the total action is the sum of quadratic action for P{z\ — z 2 ) and the fermionic action: 


S = S* + S P , 

S* = f d 2 z i, 2 $ at (;?i) 


< 5 (^i - z 2 )Y'd,i + S(zi - z 2 ) m + I Q J P(zi - z 2 )G{z l - z 2 ) 


0 0 




S P = N J d 2 zi )2 G(zi - z 2 )P(z l - z 2 )P(z 2 - zi). 

In momentum representation 

S* = J d 2 q^(-q,-q)K$ a (q,q), 

K=t m T( ) _ f d 2 , P(k,k) 

q + T(q ) m) ’ [q) J k (q - k ) 2 ’ 

" ' 2 „, 2 ] P(k,k)P(q,q) 


(A.7) 


(A.8) 


S P = N / d z / d z qd z k 


{q - ky 


Integrating over fermions gives 


Z* = J [c/T]e“ s ' 1 ' = e 


-Sv _ e N Tr log K 


(A.9) 


where 


Tr log K = / d 2 zd 2 q log det/i(g, q) — I d z zd z q log ( m 2 — |g| 2 — q I d z k 


2 „ j2. 


Combining (A.10) with Sp from (A.8) we obtain the total action for P(q,q ): 




S = N I d 2 zd 2 q 


2 u P(k,k)P(q,q) , / 2 .12 ^ , 2, p ( k , k ) 


d z k 


(q - 


— log m — |g| ~q d k 


(q - k)' : 


(q — k ) 2 

(A.10) 


(A.ll) 


In the large-A T limit solution to the theory is obtained simply by extremization of (A.ll) 
w.r.t. P(q): 


SS — N / d z zd z qd z k l 


2_j2 i2t. 


(q - h ) 2 

The corresponding equation of motion is 

P{q,q) = 


P{q,q) 


\q\ 2 - m 2- q J d 2 k 


(A-12) 


g| 2 -m2-g/d2 fc ^(M) 


(A.13) 


(* 2 ), 
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We can see that the equation (A. 13) is just’t Hooft gap equation [11] on full fermionic 
propagator P(q,q ) in leading order in large-TV limit. Indeed, the full fermionic propagator 
is equal to the sum of geometric series of 1PI propagators. On the other hand, 1PI 
propagator is equal to [11] 

S M = / d 2 k ^'*y 2 (A.14) 

and therefore the sum of geometric series is equal to 


1 | g | 2 -m 2 S (^W) 

This agrees with the equation (A. 13). 


i 2 -iS lpk j^w 


(A.15) 


Appendix B. Bethe-Salpeter equation for four-quark vertex 

Let us give a short review of how one can arrive at the Bethe-Saltpeter equation (5.31) 
on four-quark vertex T. Consider four quarks scattering process, with two incoming quarks 

r'fc)) (b.i) 

where oq is a fixed color, and two outgoing quarks 

(> 3 (p 3 ), ^ Q3 *CP4)) • (B.2) 

Denote q — pi — p 2 to be momentum flowing in this process. Due to (5.21) such diagram 
is equal to 

8{pi -P2+P3- P4)P(Pl)P(P2)P(P3)P(P4)T(p 1 ,p 3 -, q) • (B.3) 

On the other hand [12], this diagram is equal to sum of the diagram 

. 2 

S{pi ~P2 +P3 -P4)P{P1)P{P2)P{P3)P{P4)-, -—-72 , (B.4) 

[Pi ~P3) Z 

describing scattering process in which incoming quarks (B.I) exchange one gluon and fly 
away as (B.2); and the diagram 

J d 2 k ig 2 P(k)P(k-q ) - _* T(k, p 3 ; q)P(p 3 )P ( Pi ), 

(B.5) 

describing process in which two incoming quarks (B.I) exchange one gluon with momentum 
k — pi , scatter in two outgoing quarks 

$<(kj) , (B.6) 

which subsequently scatter on vertex T(k,p 3 ;q), giving two outgoing quarks (B.2). This 
way we obtain equation (5.31). 


5(pi-P2+P3~P4) P(pi)P(p2) 
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